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ABSTRACT 


A steady capillary driven flow is developed for a liquid index in a circular tube which is 
partially coated with a surface modifier to produce a discontinuous wetting condition 
from one side of the tube to the other. The bulk flow is novel in that it is truly 
steady, and controlled solely by the physics associated with dynamic wetting. The 
influence of gravity on the flow is minimized through the use of small diameter tubes 
~ O(lmm) tested horizontally in a laboratory and larger tubes ~ O(lOrara) tested in 
the low gravity environment of a drop tower. Average steady velocities are predicted 
and compared against a large experimental data set which includes the effects of tube 
dimensions and fluid properties. The sensitivity of the velocity to surface cleanliness 
is dramatic and the advantages of experimentation in a microgravity environment are 
discussed. 


I. Introduction 

The capillary driven flows to be discussed are constrained to capillary flows where 
spatial gradients in the interfacial tension, cr, are absent (V • o = 0). Such flows exist 
solely as a result of system geometry and surface wettability and though they are 
ubiquitous in nature and industry they are particularly conspicuous in low gravity 
environments such as those produced aboard orbiting space craft. In this study, a 
simple experimental technique is adopted which permits the observation of purely 
capillary driven flows which are steady. Such flows may prove useful as a tool to 
study fundamental characteristics of capillary flows relevant to many fluids manage- 
ment processes. The reported experiments are intended to be ‘demonstrative’ of this 
possibility. 

In fig. 1 a right circular cylindrical tube partially filled with a liquid index is 
sketched in a zero gravity environment. The ends of the tube are open to the atmo- 
sphere and the left hand side of the tube is treated such that a less wetting condition 
(higher contact angle, 0) exists there with respect to the right hand side (0i < 02). 
Since the pressure jump across a given interface is proportional to the curvature of 
that interface, a pressure gradient in the direction of the more favorably wetted side of 
the tube is established which, barring the effects of contact angle hysteresis, produces 
a flow resisted by the fluid’s viscosity. After a brief transient, the macroscopic flow is 
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steady and persists as long as the liquid slug bridges the wetting discontinuity. 

The earliest reference found explicitly envisioning such a flow (though there are 
likely to be earlier works) is that of Reynolds and Satterlee [1] who briefly discuss 
the case of capillary flow in tubes where the quantity cr cos 6/r is not equal at both 
interfaces, r being the local tube radius. One special case is that of Fig. lb which 
happens to produce a flow which is steady. This feature is unique in two respects: (1) 
It can avoid the need for transient analysis common, for instance, in similar “capillary 
rise” experiments [2] and (2) It is still a true capillary driven flow in which the steady 
fluid velocity is not maintained by, say, pump pressure, gravity head, a plunging tape, 
etc., but rather by the dynamic interface curvature alone. 

To demonstrate the technique, a 1-dimensional analysis is performed to determine 
the steady average velocity, V, of a liquid index of length l, viscosity /x, density p, and 
surface tension a, in a tube of radius r, with air as the host medium. Simple experi- 
ments are then outlined and data presented for tests performed using small capillary 
tubes (r < 1.0mm) in a laboratory and larger capillary tubes (r < 5.0mm) in a drop 
tower. The experimental results are scaled and presented against the independent 
variables of the problem. A brief discussion follows. It is no surprise that these flows 
are quite sensitive to the level of cleanliness of the tubes and that reproducibility and 
predictability are greatly increased with increased control of the surface conditions. 
Such sensitivity to surface conditions may make these flows ideally suited for studies 
of capillary flows over heterogeneous surfaces where the “heterogeneous” nature of 
the surface is highly controlled. For examples of such studies see [3-4]. 


II. 1-D Analysis; Average Steady Velocity 

A sketch of the slightly generalized problem of interest is provided as fig. 2. The 
principal assumption, evident in the figure, is that the menisci can be approximated 
as spherical caps. This requires that the Weber, Bond, and Capillary numbers are 
small, We = prV 2 /<r < 1, Bo r = pgr 2 /a <C 1, and Ca = Vp/cr < 1, where g is the 
acceleration field strength. A 1-dimensional analysis follows naturally incorporating 
the further assumptions of a quasi steady Poiseuille velocity profile, vtjr 2 » 1, where 
v — p/p and t is time, a slug length such that r/l 1 (reducing the need to account 
for 2-dimensional effects in the vicinity of the menisci), constant properties, slug 
length, tube radius, and the use of a dynamic contact angle simplification. With this, 
the resulting force balance is written directly 

m~ = F„ - + F g (1) 

where the forces due to surface tension, viscosity, and gravity are F a = 2irrcr(cos $i - 
cos 0 2 ), Fn = 8?r plV + 8?r p g (L - l)V, and F g = 7rpr 2 ^sin0, respectively, m is the 
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mass of the slug, 0i and 0 2 are the “apparent” dynamic advancing and receding 
contact angles, 1 and 4 > is the tube tilt angle with respect to the tube axis. fi g is 
the gas viscosity and L is the overall tube length. The horizontal is defined by the 
orientation of <7, and (j) is chosen positive in the counter clockwise direction. In the 
limit of low Capillary number the apparent dynamic contact angle is often defined in 
the literature as the geometric angle determined by extrapolating the bulk meniscus 
curvature to the wall. Angles defined as such are functions of Ca and are often 
determined by empirical techniques such as that used by Hoffman [5]. The quantity 
cos 0i — cos 0 2 = A cos 0 is thus a function of V and the steady value of V can be solved 
iteratively with empirical data for 0i and 62 as functions of Ca. A less empirically 
burdensome result for V may be obtained by assuming A cos 0 « const. This choice 
does not alter the value of V at long times and allows the integration of eq. (1) which 
gives 

V — f A cos 0 4 - sin <j> ] (2) 

4 (lI \ 2 r ) 

for the steady average velocity with 


M = 1 + — 



Entrance/exit region effects in the gas phase have been ignored. The body force and 
transient terms may be ignored by proper selection of fluid properties, tube Bond 
number, and tilt angle such that 


Bo r l sin (j> 

2 r A cos 0 ^ 


(3) 


With these further simplifications, eq. (2) reduces to 

V = 4-tM- 1 Ac os6 (4) 

4 fi l 

which is the solution to Washburn’s equation [6] where the driving pressure gradient 
is given by A P = 2crAcos0/r. 


A similar analysis applied to a U-shaped circular tube, pictured in fig. 3, can also 
be performed for a ‘U-tube’ of U-bend diameter D u (the subscript u refering to the U- 
tube configuration, refer to figure for notation). Assuming the menisci reside at z = 0 
at t — 0, the body force term of eq. (1) is modified, F g — irpr 2 g ( D u cos 4> + 2/isin 0), 
and when back-substituted into (1) and integrated over time produces 


a t 


1 ( A cos 0 4- — u cos 4 > 

\ 2 r 



gr 


16 v 2 M 2 l 


. . at Bo r . A 

Sm ^I^lW Sm(1, (5) 


1 Use of the dynamic contact angle simplification is made with reluctance and though it captures 
much of the physics manifest at the contact line its use signals the ‘demonstrative’ nature of this 

work. 
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for the velocity of the liquid slug. The effects of tube curvature on the presumed 
Poiseuille flow profile are not considered, h = h(t) is the location of the advancing 
meniscus as measured from z = 0 and the result of eq. (5) is appropriate in the limits 


gr 4 

8 v 2 M 2 l 


sin <j> <C 


1 , 


gr 2 t 

AvMl 


sin <j> <C 1, 


and 


SMvt 


>1 


The behavior of V u is not steady, as was the case for the straight tube, unless Bo r 
or (j> = 0, 2 at which point V u reduces to V of eq. (4). 


( 6 ) 

= 0 , 


For large r it is clear that the third constraint of eq. (6) will be difficult to satisfy 
for small t. When this is the case, inertia plays a significant role in the initial transient 
behavior of the flow as well. The two time scales are identified below and characterize 
the transients associated with viscous diffusion and inertia, respectively: 



An effective damping ratio can be defined by 



The time scales of eq.s (7) and (8) are not mutually exclusive since, provided l and L 
are of sufficient size, t ~ r 2 /8Mv will eventually dominate regardless of the value of 
Therefore, the system transient for the case ( 2 < 1 will consist of the combined 
effects of short term inertia giving way to longer term viscous diffusion. When C 2 1, 
inertia plays no role and the system transient is t ~ r 2 /8Mv. These scales can be 
used effectively to size potential test configurations. 


The complete form for h(t) for the U-tube configuration is determined by noting 
V u = dh/dt and by integrating eq. (1) with h( 0) = K(0) = 0 to give the dimensionless 


result 

= (a, - 02 ) sin0 ( eXpla2t ‘ ] - 1 ~ ^ (eXpM 1 " 1} ) 

(10) 

where 

J - Bo r l{ Ac ° Se+ 2 r C ° S0 ) 



2 This analysis does not (but could) treat the flow when the trailing meniscus enters the ‘bend’ 
in the U-portion of the tube. 
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and 


k — 


gr A 

v 2 MH 


with h* = h/l and t* = vMt/r 2 . Fig. 4 shows a plot of h*(t*) and V*(t*) determined 
from eq. (10) and its derivative, respectively, for the the case of a lc5 silicone oil 
index with l — 131m? n, r — 577M71, and <p — — 7r/2 (The test data for this case will 
be presented in fig. 10). Several gravity levels are given which are assumed constant. 
The qualitative nature of the solution reveals a first accelerating and then decelerating 
behavior of the bulk flow which for this parametric range is quite sensitive to gravity 


level. 


From equations (2) and (5), the capillary flow velocity resulting from a discrete 
change in surface wettability is to a large degree proportional to r/l and, thus, the 
flow rate Q is proportional tor 3 /l. If the constraining conditions on eq. (2) or (5) are 
satisfied, the r 3 dependence of the flow rate implies that substantial increases in Q 
can be realized in low-g environments where r ~ 0(1 cm) as opposed to the terrestrial 
environment where r ~ O(lmm). 


III. Experimental 

Two tube configurations were selected for experimentation purposes to maximize 
the range of tube internal diameters which could be tested. Straight capillary tubes 
were used to observe the smaller tube diameter range whereas the U-tube configu- 
ration was used to observe the larger diameters. The U-tube tests were performed 
in the 2.2 second drop tower at NASA’s Lewis Research Center. The relevant tube 
dimensions are given in table 1. The specific experimental procedures for each test 
will be described shortly. 

For every test a portion of the tube interior was coated with the surface modifier 
FC-723 to yield the variation in surface wettability necessary to drive the flow. This 
was accomplished by partially dip-coating the tubes in a bath of FC-723 and allowing 
ample time to dry. As specified by the manufacturer of the coating (3M Corporation), 
a flourochemical polymer surface coat approximately 0.25 fxm in thickness could be 
readily applied to the interior of a desired portion of the tubes. The coating dried 
rapidly and showed no measurable traces of aging beyond several hours of room 
temperature cure. In this manner an essentially discontinuous wetting condition 
was achieved at the coating boundary. 3 The uncoated portions of the tubes were 
completely wetted by the silicone oils, static contact angle 9\ — 0. 

3 A more gradual variation in surface wettability could be produced by promoting evaporation at 
the coating boundary during the coating process. In fact, an essentially continuous change in contact 
angle can be established with a similar yet more controlled process. For an interesting illustration 
of this technique see [71- 
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The static advancing and receding contact angles for the test fluids were deter- 
mined in straight capillary tubes coated entirely with FC-723. These results along 
with the relevant fluid properties of the silicone oil test fluids are provided in table 2. 
The advancing values are included to provide the hysteresis band for the fluid/solid 
pairs. The static hysteresis band could be determined by solving equation (2) for 
A cos 0 with V = 0 yielding 


BOr l . , a ~ 

— — sm© = A cos u 
2 r 

which gives the maximum liquid mass which contact angle hysteresis will support. <t> 
was increased to the point at which flow was “almost” induced. At this point Qadv 
and Orec were inferred via the geometric relationship 

T 


9 = sin -1 


cr 


1 -f cr 


with a = Zi/r , being the meniscus height at the advancing (i = 1) and receding 
{i — 2) menisci (refer fig. 2). This was also the method used for the dynamic contact 
angle measurements, out of necessity, to accommodate for interface distortions due 
to mismatched tube/fluid refractive indices. 


1. Straight Tube Tests. A given straight capillary tube, once cleaned and coated, 
was charged with a liquid slug of prescribed length introduced from the coated end 
of the tube. The tube was then tilted with respect to gravity until the slug began 
to flow towards the boundary of the coating. When this boundary was reached, the 
fluid rapidly achieved a steady velocity in the direction of the uncoated side of the 
tube. At this point the tube was immediately placed on a level glass table, the entire 
process being backlit and filmed with a variable speed cine camera directed vertically 
downward at the tube. Due to the magnification of the camera required to resolve the 
meniscus (< 8 mm FOV) and the long slug lengths tested (up to ~ 300mm), it was 
necessary to move some of the tubes during the flow from the advancing meniscus 
to the receding meniscus. This “move” was made as a discrete step between ample 
periods of cine camera photography on the advancing and receding meniscus regions. 
The glass table was leveled to the degree necessary to satisfy eq. (3). This required the 
use of a precision level and grades of less than 170 {xm/m (cj) < 0.01°) were maintained. 

The uncertainties particular to the straight tube tests and unaccounted for in the 
development of eq. (2) are the transients required to establish a fully developed flow, 
entrance/exit region effects, and non-axisymmetric menisci attributable to the gravity 
field being perpendicular to the tube axis. For the typical test fluids, tube radii, slug 
lengths, and the resulting flow velocities r*j 0(lcm/s), these effects could be estimated 
separately and shown to have little effect (< ±5%) on the steady velocities measured. 


2. U-tube Tests. The U-tube approach was adopted as a simple means to maximize 
the 2.2s of low gravity time afforded by the 23m drop tower. The initial condition 
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was such that the liquid slug strattled the coating boundary prior to the release of 
the experiment apparatus into free fall. The step reduction in gravity marked the 
onset of capillary driven flow. Fig. 5 shows a selection of frames for two such tests. 
Backlighting and high speed cine camera photography were employed for the data 
collection. 

In addition to the entrance/exit region and transient effects, mentioned previously 
for the straight tubes, three further causes for increased uncertainty were introduced 
by the U-tube approach. These were centrifugal effects, a tube eccentricity effect 
at the bend due to their fabrication, and an additional transient necessary for the 
normal gravity menisci to establish their low-g surface shapes. Each of these were 
examined independently and a cumulative error estimated to be no more than ±10% 
of V u . 


Tube Interior Surface Conditions . As has been shown by numerous investigators, 
the interior surface condition of the tubes plays a most important role in capillary 
driven flows. To illuminate such effects, a number of tube cleaning/preparation tests 
were conducted similar to those of Mumley et al[8]. These tests included initially dry 
tubes (D), tubes which were prewet by the test fluid leaving a macroscopic film on 
the uncoated side of the tube (PW), and tubes which were prewet and then allowed 
to drain dry for several days (PD)-what Mumley et al. refer to as a prewet-dried 
condition. Two types of dry tube tests (D) were conducted: the “as delivered,” or 
factory cleaned condition, where the tubes where not cleaned but tested in the condi- 
tion as received from the manufacturer, and the “cleaned,” or lab cleaned condition, 
which consisted of a series of rinse/washes in a soap solution, distilled water, acetone, 
methanol and/or ethanol, and a distilled water flush with air dry. A final cleaning 
step was attempted by baking out the cleaned tubes in a vacuum oven, but such tubes 
could not be coated easily by the FC-723 in a discontinuous manner. This may be 
in part due to increased surface diffusion of the volatile coating. Sporadic, unsteady, 
and often nonexistent flows resulted in such cases. The U-tubes were not tested using 
the PW or PD preparatory procedures. 

The larger inner diameter U-tubes permitted direct contact measurements of the 
internal surface conditions. A lab standard surface profilometer was employed for this 
purpose and, as a demonstration, values for the mean, rms, and peak-to-valley sur- 
face roughness were found to be 0.04, 0.048, and 0.203 \im, respectively, as measured 
parallel to the tube axis over a 0.8mm sampling length (all values ±50%). Circum- 
ferential measurements could also be made. Such information is rarely available with 
experimental data on flow in capillary tubes. 
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IV. Results 


The axial location of a given meniscus was first tracked in time using a motion 
analyzer. Fig. 6 illustrates some results for a number of factory cleaned straight tube 
tests performed with lcS silicone oil and a variety of slug lengths. As seen from the 
figure, the base flow is indeed steady. Time t = 0 here represents the initiation of 
data collection and not the beginning of the flow. 

Deviation plots between instantaneous velocities measured experimentally and 
those determined by a linear fit to the data show agreement to within 1.5% or better, 
an agreement which is at or below the resolution accuracy of the framing rate of 
the cine camera. Fig. 7 provides two examples. The raw data sets were acquired 
using an automated image analysis technique which digitized the cine film records 
frame by frame and tracked the menisci across the field of view. The grabbed images 
were low-pass filtered and a median threshold was selected to define the interface. 
Only data points are shown in the figure (no fine). The linear fits to the data are 
subtracted from h measured experimentally and the differences are shown. In fig. 
7a, a worst case, the long period oscillation reveals a peak uncertainty of 1.5% for 
the instantaneous velocity and may reflect a slightly varying film speed. In fig. 7b 
no trend of the deviations is visible and discrepancies are far less than ±1%. The 
general form of a given deviation plot was neither reproducible, nor traceable to a 
given surface condition or test configuration. 

Fig. 8 displays sample U-tube data acquired similarly for several test cases where 
r 2 /SMi v < 0.25s. The initial transient is readily seen and the flows can be said to 
be predominately steady. The fact that a “falling off” of the velocity is present for 
these tests indicates the significance of the body force term, atBo r sin (j)/ A\ jlMI in 
eq. (5). Given that the low-g environment of the drop tower is ~ O(l0~ 4 g o ) due to 
aerodynamic drag ( g Q = 9.8 m/s 2 ), taking <f> = — tt/2, the principal direction of the 
drag force, the above term is the appropriate magnitude (~ 0.05t) to produce this 
behavior. This effect was also observed to diminish with decreasing time, decreasing 
tube radius, increasing fluid viscosity, and increasing slug length. These trends are in 
accord with eq. (5). Fig. 9 plots the instantaneous velocity against time for the data 
of fig. 8. The negative linear trend is obvious. The prefactor G — r 2 / 1 4 is introduced 
on the figure as a measure of the intensity of the transient term in eq. (5), where 
aBo r j\Ll — Gg. It is apparent that the slope of the linear fits to the data are related 
to G as anticipated by eq. (5). For fixed g one could minimize gravitational effects for 
the U-tube configuration by minimizing G. It should be emphasized that the straight 
tube configuration adopted for use in the drop tower would eliminate this falling-off 
effect and that tubes with L and l of sufficient length would achieve steady flows of 
similar calibre to the straight tube tests of this study. 
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Several of the larger inner diameter U-tube tests were performed where r 2 /8vM ~ 
3 s, which implied viscous time scales larger than the drop time available. An example 
is provided in fig. 10 for 1 cS silicone oil in a 9.88mm inner diameter tube. The 
initial oscillation of the interface due to the step reduction in g is more exaggerated 
than for the cases of fig. 8. The decay time for these oscillations (< 0.25s) is 
half that predicted by a recent correlation developed for a similar step reduction in 
g for partially fluid-filled cylindrical containers [9]. The discrepancy is due to the 
additional damping influence of the bulk flow. Unert for this case is approaching Is 
and C 2 — 0.112. The basic elements of the greater unsteadiness of the flow can be 
gleaned from eq. (10) which reveals a first accelerating and then decelerating flow 
due to the presence of finite g with 4> — — 7r/2 (refer to qualitative result of fig. 
4). However, the analysis cannot capture the quantitative transients of these flows. 
This should be expected since A cos# = const, is hardly representative of the true 
phenomena, the constant low -g assumption for the drop tower tests is not correct, 
and the constraint 8Mvtjr 2 ^ 1 is not satisfied. 

The distance-time data were reduced to give average velocities using a central 
difference scheme on data such as those presented in fig. 6 and 8. At least 15 values 
were used in each average value. An oddity of some of the experiments was that 
on occasion the average velocity values determined at the advancing and receding 
menisci for the same liquid slug were sometimes detected to differ by as much as 
±5%. This was never true for the U-tube tests and no pattern, or test condition or 
configuration could be established to determine when or why this would happen. In 
the event of such an occurrence the advancing and receding velocities were simply 
averaged to determine V. 

Dynamic contact angle measurements were made to within ±3° for the wide range 
of tests performed, and in many cases ±2° or better was common, the difference com- 
ing from a trade-off between accurate contact angle measurement and accurate ve- 
locity measurement. Though the point becomes mute at low Ca where measurement 
uncertainties of V were less than ±1%, the higher uncertainty of ±3° for 6 incurred 
errors up to ±25% in the calculation of V from eq. (4) through A cos 9. Cumulative 
uncertainties as high as ±30% for V and V u calculated from the measured quantities 
could be expected in these cases. Of the many cases where contact angles could be 
determined within ±2°, uncertainties could be banded between ±25%. 


Eq. (4) can be rearranged to give 


V 

A cos 0 


= — 


4/x l 


= F 


(ii) 


which separates the dependent from the independent parameters of the problem. The 
former is plotted against the latter in fig. 11 on a log-log scale. This figure represents 
the sum total of the study for all of the tests conducted. Eq. (11) is displayed with 
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an error band allowing for the estimated 25% uncertainty. Fig. 12 gives a breakout 
of the different tube surface conditions. As can be viewed in fig. 12a, and as will be 
seen shortly, not much control was afforded by the added effort of further cleaning 
the tubes from the state in which they were received. The cleaning for the U-tube 
tests (fig. 12b) was unspecified and consisted of an unknown combination of factory 
and lab clean conditions. The pre-existing films in the prewet tests (fig. 12c) for the 
straight tubes were not closely controlled or measured. By comparison, the case of 
the prewet and dried tubes (fig 12d) was highly repeatable and the data generally fell 
within the error band of eq. (11). 

The standard deviation for the repeatability of the measured quantity V/A cos 9 
for the straight tube tests was ±41% for the factory cleaned tubes, ±38% for the lab 
cleaned tubes, and ±4% for the PD tubes with F held approximately constant. Fig. 
13 presents straight tube data indicating the tube cleaning technique for each data 
point. As seen from the figure, the cleaning technique clearly influences primarily the 
degree of scatter of the data through A cos 9 , the PD condition exhibiting the most 
control. Too little U-tube test data was collected to permit reasonable statistical 
analysis. 

Rearranging eq. (11) further gives 

Ca = ^AT 1 A cos 9 = F 2 (12) 

41 

which is presented in fig. 14 using the data of fig. 11. Worthy of mention is that 
for purely capillary driven flows as defined herein, inspection of eq. (12) shows that 
the Capillary number scaled on the correct velocity, V, is a function of the geometry 
of the system (r/l) and the dynamic surface wettability (A cos 9). It should also 
be pointed out that the agreement of the data with eq. (11) is increased with the 
inclusion of M which accounts for the viscosity of the gas phase. No effort was made 
to vary M specifically and values between 1 and 1.91 were achieved. 

Though the flow is relatively well represented by eq. (11), only modest correlation 
between 9adv and Ca was established. Fig. 15 is provided as an illustrative example 
for the factory cleaned straight tube tests. Symbols are provided indicating tube inner 
diameter (fig. 15a) and fluid viscosity (fig. 15b). The correlation given by Hoffman 
[5] for Ca < 0.1 for the advancing meniscus 

9adv = 4.54 Ca 0353 (radians) (13) 

is provided for reference. Without this curve sketched on the figures the global trend 
for 9adv identified by Hoffman might not be distinguished. This is in part due to the 
limited range of Capillary numbers accessible by purely capillary driven flows. The 
data for the PD tests show little improvement. The receding contact angle values as 
functions of Ca are included as fig. 15c for completeness. 
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Fig. 16 provides the results for A cos# with Ca for the data of fig. 15 which 
shows that, in this study, the larger values of A cos 9 were achieved for the smaller 
tube radii and higher viscosity fluids. 

V. Discussion 

The scatter of the data displayed in fig. 11, and more specifically in fig. 13, shows 
the effects of surface conditions on the value of the steady capillary velocity of liquids 
in tubes with wetting discontinuities. It must be bom in mind that in most cases each 
average velocity measurement was made to a relatively high degree of precision, often 
within 1% or better. Thus, as stated, the principal uncertainty is mainly attributed to 
the contact angle measurement uncertainty of ±2° leading to an uncertainty in A cos 9 
of up to ±25%. Albeit a large uncertainty, the experimentally determined spread of 
V/A cos 0 for F fixed is even larger, up to ±41%, for the cases of the factory and 
laboratory cleaned tubes. Nonetheless, the spread is still somewhat centered about 
eq. (11) implying that the physical effects of surface conditions of microscopic scale at 
the contact line are largely accommodated by changes in interface curvature observed 
on a macroscopic scale. 

Fig. 17 displays the dynamic contact angle data as a function of Ca for the 
data of fig. 13. The factory and lab cleaned tubes yielded slower flows (up to three 
times slower than the PW tests) and hence lower values for A cos 9 with increased 
uncertainty. The fact that the dynamic receding contact angle actually increases with 
increasing Ca , as shown in fig. 17, is perhaps misleading, contrary to intuition, and 
deserves attention. From the figure it is apparent that the Hoffman curve gives an 
acceptable representation for the advancing contact angle. Substituting this form into 
eq. (5) and solving for 92 gives 

cos 02 = cos (4.54 Ca 0 - 353 ) - jCa (14) 

which shows 9 2 increasing with increasing Ca as is indicated in fig. 17. A similar 
expression could be obtained regardless of the static advancing and receding contact 
angles so long as an analytic form is available for 9\. Since #2 cannot increase past the 
static receding value, an upper bound on Ca may be established for such capillary 
driven flows when Al/r is known. As an example, the smallest value for 4 l Jr of this 
study was 22.5. For a static receding angle of 50°, eq. (14) when solved for Ca gives 
the peak value Ca = 0.0024, beyond which no flow could reach. This limit is indicated 
in fig. 15c as is eq. (14) for the maximum and minimum values of Al/r used in these 
tests. The data are reasonably well bounded between the predicted limits and the 
value for Ca < 0.0024 is indeed the limit for the data set (refer also to fig. 14). 

The increase of 9 2 with Ca is a result of the coupling between V and A cos 9. 
Since Ca is determined naturally by the system it is not directly controlled as is often 
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the case for studies which drive the contact line and thereby force large values of Ca 
beyond that which the flow could attain naturally by capillarity alone. The fact that 
for these flows Ca oc A cos0, it is no surprise that Ca is maximized for 0 2 —> 

A cos 0 can be scaled by its maximum value A cos 0 re / = 1 —cos 0 rec , where here 0 rec 
is the static receding contact angle and the advancing static angle has been set to 0. 
The contact angle values listed in table 2 are used to further reduce the data of fig. 16b 
and are replotted in fig. 18. For the fluid/solid/coating combinations selected for this 
study A cos 0/A cos 0 re / centered around 0.51 with a standard deviation of 35% for all 
of the tests performed. This may provide an order of magnitude predictive capability 
to estimate V for a given test configuration, given that information concerning the 
static contact angles on the test surfaces is known. 


Concluding Remarks 

The cursory experimental results presented herein thoroughly demonstrate the 
steady nature of flows generated by discontinuously wetted tubes. This is most clearly 
seen from the individual tests performed (as an example, refer to fig. 6). The general 
agreement of eq. (11) with the complete data set in fig. 11 is also indicative of the 
accuracy of the simplified analysis over a wide parametric range, the scatter of the 
data about the predicted curve being due primarily to a large uncertainty ~ ±25% 
in the measurement of A cos 0. Provided empirical data or analytic expressions are 
available for the advancing and receding contact angles as functions of Ca, eq. (4) may 
be solved directly for V. What is interesting here is the degree to which these flows are 
sensitive to tube surface conditions. For example, fig. 17 displays a three-fold increase 
in Ca for the PD condition. These results indicate the role of microscopic effects in 
the near vicinity of the contact line, since all else, including the static advancing and 
receding contact angles, is held constant. The technique of using surface coatings to 
produce steady capillary driven flows may serve well as a tool to probe the innerscale 
mechanisms controlling dynamic wetting. 

The advantage of experimentation in the low gravity environment is that ‘capil- 
lary’ tubes can be quite large. Analytic models of such effects as surface roughness, 
surface cleanliness, patchy wettability, etc. on capillary driven flows could be further 
developed utilizing this test configuration since the internal conditions of the tubes 
can be quantified in detail. Tubes could be selected large enough that controlled 
surface roughness could be accurately machined, micro-fabricated through chemical 
deposition, selectively etched, etc. Thus, the low gravity environment presents further 
possibilities which, if exploited, could further our understanding of the complicated 
nature of capillary driven flows on non-ideaiized surfaces. 
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Table 1. Tube data. 


Tube TVpe 

ID (±2%) 
(mm) 

L (±2%) 
{mm) 

D u ±X 
(mm) 

range of l tested 
(mm) 

Straight 

0.516 

95.0 

- 

6-79 

(thin-walled) 

0.730 

95.0 

- 

5-79 


1.145 

95.0 

- 

39-42 


2.312 

95.0 

- 

19-44 

Straight 


212-481 

- 

60-303 

(thick-walled) 

mmm 

102 

- 

15-20 

U-tubes 

2.29 

94, 121, 173, 206 

10.3 (0.5) 



3.81 

106, 133, 186 

17.2 (1.0) 

66-114 


5.94 

113, 192, 140 

22.0 (0.9) 

79-121 


7.77 

121, 139, 190, 206 

22.6 (0.9) 

85-133 


9.88 

158, 206 

32.2 (1.2) 

119-134 


Table 2. Properties of silicone oil test fluid on FC-723 coated glass. 


v ±2% 
(xl0 -6 m 2 /s) 

p± 5 
kg/m 3 ) 

B 

0ad v ±X 

(deg.) 

Ortc 

(deg.) 

0.65 

760 

0.0159 

44.7 (5) 

35.0 (5) 

1.0 

816 

0.0174 

40.9 (2) 

31.0 (2) 

1.5 

851 


51.7 (2) 

41.7 (2) 

2.0 

872 

0.0187 

52.2 (2) 

47.3 (2) 

5.0 

913 

0.0197 

57.4 (2) 

48.9 (1) 

10.0 

935 


59.3 (3) 

51.3 (2) 

20.0 

949 


56.0 (3) 

46.2 (2) 

50.0 

960 

0.0208 

59.3 (3) 

54.5 (2) 

100.0 

964 

0.0209 

- 

- 
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Figure 1 . — Schematic of steady flow resulting from a dicon- Figure 2. — Schematic of liquid slug of length I in a straight 

tinuous wetting condition on the interior of a circular tube capillary tube of length L. 

of constant radius where 0-j < 0 2 . 



Figure 3. — Schematic of liquid slug in a U-tube where each 
leg of the 'Ll 1 has a different wettability with the fluid. 
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Figure 4a. — h* ( and 4b. , versus t*, respectively, for 

the U-tube configuration: effect of finite acceleration 
field, g (m/s 2 ). 


Figure 5. — Selected frames from two U-tube 
drop tower tests, v = IcS. Time progression is 
top to bottom in 0.5s increments, U-portion of 
tubes not pictured, and left legs of U-tubes 
are coated. [Left: r = 2.97, 1 = 1 21 , L = 1 92 mm, 
V u = 8.1 mm/s. Right: r = 1 .90, 1 = 115, L = 

186 mm, V u = 3.2 mm/s.] 








14 


G 


e 20.3 
• 12.8 
v 15.2 
▼ 4.0 



o 0.50 1 — 


• 1.14 64 1.0 

• 1.14 101 1.0 

v 2.97 116 5.0 

▼ 1.14 65 5.0 


0.5 1.0 1 .5 2.0 2 

Time, s 

Figure 8. — Meniscus axial location versus time for 
several U-tube tests. (Lines are linear fits to the 
data above 0.25s.) 


Figure 9. — Instantaneous velocity of liquid slug versus 
time for U-tube tests of fig. 8. Symbols are consistant 
with test conditions. 
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Figure 10. — Meniscus location for U-tube test of 
large radius, v = IcS, and r = 4.94, 1 = 131 , and 
L = 206 mm. 
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Figure 1 1 . — V/Acos0 versus F for tests performed. 
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Figure 13. — Effect of tube cleaning condition for initially 
'dry' tubes. 
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Figure 15— Ca dependence of (a) 0 a dv ID indicated, 
(b) 6 a dv viscosity indicated, and (c) 0 rec viscosity 
indicated for factory cleaned (D) straight tubes. 
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Figure 17. — Dynamic contact angle versus Ca for fixed F, 
data of fig. 1 3. Hollow symbols are advancing values, 
solids are receding. (41/r = 438.) 



Figure 16. — AcosG versus Ca for factory cleaned (D) 
straight tubes, (a) ID indicated, (b) viscosity indicated. 



Ca 

Figure 18. — Reduced AcosG versus Ca for data of fig. 16b. 
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